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The process e + e~ — > 7r + 7r~7r°7 has been studied at a center-of-mass energy near the T(4S) 
resonance using a 89.3 fb -1 data sample collected with the BABAR detector at the PEP-II collider. 
From the measured 3ir mass spectrum we have obtained the products of branching fractions for the 
lj and <f> mesons, B{u -> e + e")B(uj -> 3tt) = (6.70 ± 0.06 ± 0.27) x 10~ 5 and B((f> -> e + e -)B{4> -> 
3tt) = (4.30 ± 0.08 ± 0.21) x 10~ 5 , and evaluated the e+e" -> tt+tt^tt cross section for the e + e" 
center-of-mass energy range 1.05 to 3.00 GeV. About 900 e + e~ — > J/ipj — > tv + tt~tv°j events have 
been selected and the branching fraction B(J/tp — * 7r + 7r~7r°) = (2.18 ± 0.19)% has been measured. 
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FIG. 1: The diagram for the e + e — > 7r + 7r 7r°7 process. 



INTRODUCTION 



The analysis of the process e H 



7" 



- hadrons, where 



the photon emission is caused by initial state radiation 
(ISR), can be used to measure the e + e~ annihilation 
cross section into hadrons over a wide range of center- 
of-mass (cm.) energy in a single experiment [1-3]. In 
these events, the invariant mass, ys', of the final state 
hadronic system corresponds to the "effective" center-of- 
mass energy after ISR. 

This method, currently used both at the KLOE exper- 
iment at DAFNE (Frascati) [4] and at the BABAR exper- 
iment at the PEP-II B Factory (SLAC), is applied here 
to study the process e + e~ — > tt + -k^-k° at low energies (1 
to 3 GcV). 

The Born cross section for the e + e~ — > hadrons+7 pro- 
cess (Fig. 1) integrated over the momenta of the hadrons 
is given by 



da(s, x, 9) 
dx d cos 9 



W(s,x,6)a (s(l-x)), 



(1) 



where -y/s is the e + e~ center-of-mass energy, x = 
2Ej/ -y/s, E 1 and 9 are the photon energy and polar angle 
in the cm. frame, and s(l — x) = s', already mentioned 
above. Here <tq is defined as the Born cross section for 
e + e~ — > hadrons. The so-called radiator function (see, 
for example, Ref. [3]) 



W(s,x,0) 



2x 



x 
~2 



(2) 



describes the probability of ISR photon emission for 
9 3> m e j \fs. Here a is the fine structure constant and m e 
is the electron mass. The ISR photons are emitted pre- 
dominantly at small angles relative to the initial electron 
or positron directions; however about 10% of the photons 
have cm. polar angles in the range 30° < 9 < 150°. In 



the present analysis, we require that the ISR photon is 
detected. 

The differential cross section for ISR production of a 
narrow resonance (vector meson V), such as J/ip, decay- 
ing into the final state / can be calculated using [3] 



dajs, 9) _ 12n 2 T(V ^ e+e-)B(y -*J) 
d cos 9 ms 



W{s,x ,9), 



(3) 

where m and T(V — ► e + e ) are the mass and electronic 
width of the vector meson V , xo = 1 — m 2 /s, and B(V — > 
/) is the branching fraction of V into the final state /. 
Therefore, the measurement of the number of J ftp — ► 3tt 
decays in e + e~ — > Sttj determines the product of the 
electronic width and the branching fraction: T{J /ip — > 
e+e-)B(J/i/j -► 3tt). 

The e + e~ — > 7r + 7r~7r° cross section in the energy region 
Vs 7 < 1 GeV is dominated by the w(782) and 0(1020) 
mesons 1 . This energy region has been studied in many 
experiments with high statistics and the Particle Data 
Group (PDG) parameters [5] for w and <f> mesons have 



relatively high precision (2-3% for B(V 



-)B(V 



37r) and about 1% for the total widths). 

The energy region above the <fi was studied in two ex- 
periments: SND [6] for energies up to 1.4 GeV with sta- 
tistical precision about 10% and DM2 [7] for energies in 
the 1.34-2.40 GeV range with statistical precision about 
25%. As pointed out in Ref. [6], there is a significant sys- 
tematic shift between these two datasets, and the DM2 
data need to be scaled by a factor 1.72 ± 0.24 in order to 
fit with those of SND. 

In this energy region, the e + e~ — * 3-7T cross section is 
generally described as the sum of two resonances u/(1420) 
and £j"(1650). So cross section measurement allows the 
determination of the uj' and u" parameters. Masses, 
widths, and decay modes for these resonances are not 
well established. The PDG [5] gives only estimates for 
these parameters. 

The main goal of this analysis is an independent mea- 
surement of the e + e~ — > 7r + 7r~7r° cross section in the 
energy region from 1.05 to 3.00 GcV. The aim is to sig- 
nificantly improve the precision of the cross section for 
energies above 1.4 GeV. Our data in the u> — <f> region 
can be compared with the more precise e + e _ data in 
this region and the difference can be used to check our 
systematic error estimation. 

We study J/ip production in the process e + e~ — > 
7r + 7r~7r°7, and measure the product T(J/ip — > 
e + e~)B(J/'4' — > 3?r). The branching fraction B(J/tp — > 
3n) is then determined using the known value of 
T(J/tp — > e + e~) [8]. The decay J/ip — > 37r has been stud- 
ied in many experiments [9-15] but only three of them 
measured the decay rate without any restrictions on the 
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Throughout this paper, 2tt, 3ir, and 4tt mean tt+tt , tt+tt 7r , 
and tt+it - TT°iT , respectively. We also use the notations u>, <j>, 
w', and w" for w(782), 0(1020), w(1420), and oj(1650). 
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invariant mass of the two-pion system. Three results for 
B(J/ip -> 3tt) are (1.6 ± 0.4)% [12], (1.5 ± 0.2)% [13], 
(1.42 ± 0.19)% [14]. The average of these measurements 
is (1.47±0.13)%. This value is in significant disagreement 
with more recent result from the BES Collaboration [16]: 
(2.10 ±0.12)%. 



II. THE BABAR DETECTOR AND DATA 
SAMPLES 

In this paper a data sample of 89.3 fb _1 , collected 
by the BABAR detector [17] at the PEP-II asymmetric- 
energy storage ring, is analyzed. At PEP-II, 9-GeV elec- 
trons collide with 3.1-GcV positrons at a center-of-mass 
energy of 10.6 GeV (T(4S) resonance). 

Charged-particle tracking for the BABAR detector is 
provided by a five-layer silicon vertex tracker (SVT) and 
a 40-layer drift chamber (DCH), operating in a 1.5-T 
axial magnetic field. The transverse momentum reso- 
lution is 0.47% at 1 GeV/c. Energies of photons and 
electrons are measured with a CsI(Tl) electromagnetic 
calorimeter (EMC) with a resolution of 3% at 1 GeV. 
Charged-particle identification is provided by ionization 
measurements in the SVT and DCH, and by an inter- 
nally reflecting ring- imaging Cherenkov detector (DIRC). 
Muons are identified in the solenoid's instrumented flux 
return, which consists of iron plates interleaved with re- 
sistive plate chambers. 

Signal and background ISR processes are simulated us- 
ing Monte Carlo (MC) event generators based on the 
computer code described in Ref. [18]. The event gener- 
ator for the e + e~ — > Stt^ reaction uses a model with an 
intermediate pir state. The background process e + e~ — > 
7r + 7r~7r 7r°7 is simulated with W7r° and ai(1260)7r inter- 
mediate states in the proportion that matches existing 
experimental data. The extra soft-photon radiation is 
generated with the use of the structure function method 
of Ref. [19] and the PHOTOS package [20] for electron 
and charged hadron brcmsstrahlung, respectively. Since 
the polar-angle distribution of the ISR photon is peaked 
near 0° and 180°, the events are generated with the re- 
striction on the photon polar angle in the cm. frame, 
20° < 6 < 160°. We also require that the invariant mass 
of the hadron system and ISR photon together is greater 
than 8 GeV/c 2 . This second cut restricts the maximum 
energy of extra photons emitted by the initial particles. 
The background processes e + e~ — ► 7r + 7r~7, /i + /^~7 are 
generated with the Phokhara program [21], which in- 
cludes next-to-leading-order QED corrections and sim- 
ulates the emission of two hard photons at large angle by 
the initial particles. The background from e + e~ — > qq 
and e+e~ -> t + t~ is simulated with JETSET [22] and 
KORALB [23] packages, respectively. The interaction of 
the generated particles with the BABAR detector and the 
detector response are simulated using the GEANT4 [24] 
package. The simulation takes into account the variation 
of the detector and accelerator conditions, and in partic- 
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FIG. 2: The x 2 distributions for data (points with error bars) 
and simulated (histogram) events from the ui mass region. 
The shaded histogram shows the distribution for simulated 
background events. 

ular describes the beam-induced background, which leads 
to the appearance of photons and tracks overlapping on 
the events of interest. 



III. EVENT SELECTION 

The initial selection of e + e~ — > ir + ir~ir°j candidates 
requires that all the final particles are detected inside 
a fiducial volume. (All kinematic variables used in the 
paper are defined in the laboratory frame unless other- 
wise stated.) Since a significant fraction of the events 
contain beam-generated spurious tracks and photons, we 
select events with two or three tracks and at least three 
photons that have energies above 100 MeV and polar an- 
gles in the range 23° < 9 < 137.5° (the corresponding 
angular range in the cm. frame is 38° < 6 < 154°). 
One of the photons is required to have an energy in 
the cm. frame above 3 GeV. Two of the tracks must 
originate from the interaction point, have a transverse 
momentum above 100 MeV/c and be in the polar angle 
region between 23° and 140°. Background events from 
the process e + e~ — * e + e~7 are suppressed by requir- 
ing the ratio of the calorimeter-deposited energy to the 
track momentum, Eemc/p, to be below 0.9 for the two 
highest-momentum tracks. 

The photon with greatest cm. energy is assumed to 
be the ISR photon. The remaining photons are paired to 
form candidate 7r°s, requiring that their invariant mass 
must be in the range 0.07 to 0.20 GeV/c 2 . A kinematic 
fit is applied to the selected event, imposing energy and 
momentum conservation, and constraining the candidate 
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FIG. 3: The \ 2 distributions for data from the mass range 
1.05 < M 37r < 3.00 GeV/c 2 . The shaded histo gram shows 
events rejected by background suppression cuts. 



7T° invariant mass. The MC simulation docs not accu- 
rately reproduce the shape of the resolution function for 
the photon energy. This leads to a difference in the dis- 
tributions of the x 2 °f the kinematic fit for data and 
simulation. To reduce this difference only the measured 
direction of the 1SR photon is used in the fit; its energy 
is a free fit parameter. In the case of events with three 
tracks, the fit uses the parameters of the two tracks with 
opposite charge that have the minimum distance from 
the interaction point in the azimuthal plane. For events 
with more than three photons all possible combinations 
of photons are tested and the one with minimum \ 2 is 
used. The events with very high x 2 (> 10000) are con- 
sidered as not reconstructed. 

The x 2 °f the kinematic fit is used to discriminate real 
3-7T7 events from background. Figure 2 shows the x 2 dis- 
tribution for events from the 3tt mass region near the 
lj mass (0.75-0.82 GeV/c 2 ), where the contribution of 
background processes is small. Events with x 2 < 40 are 
selected to analyze the 3tt mass spectrum. The rest of the 
sample (10000 > x 2 > 40) is used to study both back- 
ground processes and possible selection-efficiency correc- 
tions due to data-MC simulation differences in the x 2 dis- 
tribution. The x 2 distribution for masses 1.05 < < 
3.00 GeV/c 2 (Fig. 3) shows that a significant fraction 
of events in this mass range correspond to background 
processes. 

The main sources of background are other ISR pro- 



cesses (e 



7T + 7T 7r°7r°7, 



7T+7T 7, K + K 7r°7, etc.) 



annihilation to qq and t + t~. Additional 
background suppression cuts are applied to improve the 
signal-to-background ratio in the mass region of interest. 
The ISR events with kaons in the final state (e + e~ — > 



7r" energy. 

e + e~ — ► 7r + 7r~7 and e + e~ 



K + K~it°j, e + e~ — > K + K~j) arc suppressed using the 
kaon identification based on dE/dx measurements in the 
tracking devices, and the value of the Cherenkov angle 
and the number of photons measured in the DIRC. The 
requirement that none of the charged tracks is identified 
as a kaon rejects about 95% of the kaon-induced back- 
ground with only 4% loss of 3-7T7 events. 

The radiative events e + e~ — > 7r + 7r~7 and e + e~ — > 
p + p~j with extra photons having 77 invariant mass 
close to that of a ir° arc suppressed by a cut on the 
The cut E^a > 0.4 GeV rejects 80% of 
p + p~j events with about 
4% loss of 37T7 events. 

The process e + e~ — > 7r + 7r~7r 7r°7 with a soft tt° is the 
main source of background for the process under study. 
Some fraction of 47r7 events reconstructed as 7r + 7r _ 7r°7 
contain a tt° among extra photons. We select these events 
by performing a kinematic fit for the 47r7 hypothesis. The 
requirement xln-y > 40 rejects about 40% of 47r7 events 
and only 2% of 3-7T7 events. 

The main source of background from e + e~ — > t + t~ 
is the events in which both r's decay into 27tza The 
hard photon arises from a ir° decay. Since the r — + 2irv 
decay proceeds mainly via pis, such events must have 
the invariant mass of the most energetic photon and one 
of the charged pions peaked near the p mass. The cut 
M-Tr-y > 1.5 GeV/c 2 almost fully rejects e + e~ — ► r + r~ 
background and leads to only 0.3% loss of 37r7 events. 
The remaining t + t~ background is estimated to be less 
than 0.1% of 3777 events. 

Another possible background source is events from 
e + e~ annihilation into hadrons containing a very ener- 
getic 7T°. A fraction of these events can be seen in the 
distribution of invariant mass (M* ) of two photons, one 
of which is the most energetic in an event. The second 
photon in the pair is required to have an energy above 
100 MeV. Once all possible photon pair combinations 
are checked, the one with closest invariant mass to the 
ir° mass is chosen. The events with 0.10 < M* < 0.17 
GeV/c 2 are rejected. 

The x 2 distribution for all rejected events is shown as 
the shaded histogram in Fig. 3. It is seen from Fig. 3 that 
we reject more than 60% of background events without 
significant loss of signal events. 



IV. BACKGROUND CALCULATION AND 
SUBTRACTION 

The x 2 distribution and 3ir mass spectrum for data 

> 7r + 7r~7r°7 and background 



and for simulation of e 
processes after imposing the background suppression cuts 
are shown in Figs. 4 and 5. The remaining background 
can be divided into two classes with different x 2 distri- 
butions. The first class includes e + e _ — > K + K~tt°^ and 
e + e~ — » 7r + 7r~7r 7r processes, which have x 2 distribu- 
tions peaked at low x 2 . The second class includes all 
other background processes. The simulated x 2 distribu- 
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FIG. 4: The x 2 distribution for events from the mass range 
1.05 < M3 n < 3.00 GeV/c 2 after background suppression 
cuts. The points with error bars show the data distribu- 
tion. The histogram is the sum of simulated distributions 
for e + e~ — > 7r + 7r~7r°7 and background processes. The dark 
and lightly shaded histograms show the distributions for 
e + e~ — > 7v + tt~ 7v°tt° and other background processes, respec- 
tively. 



tion for these processes is shown by the lightly shaded 
histogram in Fig. 4. The background events from the 
first class must be subtracted bin-by-bin from the 
spectrum. 

The mass distribution for kaon events surviving the 
selection (Nqk) can be obtained from the distribution of 
events with two identified kaons (N 2 k)'- Nqk{Mz-k) = 
N 2 K{M 3 T r )RK- The coefficient Rk is determined from 
the simulation, which uses kaon identification efficiency 
corrections obtained from a pure kaon sample from D* 
decays to improve the agreement of data and simulation. 
The value of Rk is found to be independent of and 
equal to 0.09 ± 0.01. The total fraction of kaon events in 
the final event sample is estimated to be 0.4%. 

The simulation shows that about 80% of the e + e~ — > 
qq events that pass the selection criteria have 7r + 7r~7r 7r° 
as the final state. In order to cross-check the value of the 
yield given by JETSET fragmentation for this particular 
final state, we use the following procedure to extract the 
mass distribution for e + e~ — * 7r + 7r~7r 7r° from experi- 
mental data. We select events with two charged particles 
and four photons with energy more than 0.1 GeV, at least 
one of them with cm. energy more than 3 GeV, perform 
a kinematic fit to the e + e~ — > 7r + 7r~7r°7r° hypothesis and 
require the \ 2 of this fit to be less than 20. The number 
of selected 47r events is found to be about 15% less than 
the number expected from JETSET. We also studied var- 
ious two- and three-particle mass distributions and find 
that both in data and simulation the e + e~ —> Tr+Tr - ^ ?? 



FIG. 5: The 3tt mass spectrum for data events with \ < 
40. The lightly shaded histogram shows the mass spectrum 
for calculated background. The dark histogram is e + e" — > 



7r + 7r 7r°7r° background. 



process proceeds via the p7T7r state. No other inter- 
mediate states are seen. To obtain the 37r mass spec- 
trum for e + e~ — * tt + it~ tt°tt events reconstructed in 
the e + e~ — ► 7r + 7r~7r°7 hypothesis ((dN/dm)^^), we 
multiply the 37r mass spectrum of selected 47r events 
((dN / dm)^) by the ratio of corresponding simulated 



distributions: (dN/dm)^ / (dN/dm)^' . The resulting 
mass spectrum is shown in Fig. 5 as the dark shaded 
histogram. The fraction of this background does not ex- 
ceed 10% in the mass region below 2 GeV/c 2 and rises 
at higher masses. 

The main contribution to the second class of back- 
ground events comes from the e + e~ — > 7r + 7r~7r°7r°7 and 
e + e~ — > 7r + 7r~7 processes. To estimate the accuracy of 
the MC simulation prediction of the background level for 
these processes we use events with x 2 > 40. Fitting dis- 
tributions of the 37T and 7r + 7r~ invariant masses for these 
predominantly background events with a sum of distri- 
butions for the signal and background processes, we find 
scale factors of (1.00 ± 0.25) for 47r 7 and (3 ± 1) for 27r 7 . 
Figures 6 and 7 show the fitted and M^+^- distri- 
butions for events with 100 < x 2 < 500. The quoted 
errors in the scale factors are much larger than the sta- 
tistical errors from the fits and originate from systematic 
effects. For 47r7 events the error is determined by the de- 
pendence of the scale factor value on % 2 . For 2ir~f events 
we observe a significant difference in the shape of the 
37r mass distribution between data and simulation and 
the 30% systematic error has been assigned to cover the 
uncertainty on the mass spectrum. 

We attribute this difference and the large value of the 
27T7 scale factor to the inaccuracy of the simulation of 
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FIG. 6: The M 3w spectrum for events with 100 < \ < 500 
and E n o > 1 GeV. The points with error bars are data. 
The histogram is the sum of spectra for simulated events of 
e + e~ — > 7r + 7r _ 7r°7 and background processes. The lightly 
shaded histogram is the sum of all background processes. 
The difference between the dark and lightly shaded his- 
tograms shows the e + e~ — > 7r + 7r~7r 7r°7 contribution. 



FIG. 7: n ir invariant mass spectrum for events with 
100 < x 2 < 500 and 0.85 < M 3v < 1.05 GeV/c 2 . The 
points with error bars are data. The histogram is the sum of 
spectra for simulated events of e + e~ — > 7r + 7r _ 7r°7 and back- 
ground processes. The lightly shaded histogram is the sum 
of all background processes. The difference between the dark 
and lightly shaded histogram shows the e + e~ — > 7r + 7r~7 
contribution. 



7r- meson nuclear interactions in the EMC. This leads to 
a difference between data and MC simulation in the mul- 
tiplicities of spurious photons arising from the nuclear 
interaction, and their energy and angular distributions. 
This conclusion is supported by the fact that the scale 
factor for e + e~ — ► /i + /i~7 background events is found to 
be close to unity. The e + e~ — ► /i + /i~7 events were se- 
lected using muon identification criteria based on EMC 
and IFR information. The level of e + e~ — > /i + /i~7 back- 
ground is found to be (5-15)% of the 27r7 background. 

Another source of the second class of background, non- 
4.7T e + e~ — > qq, q = u, d, s, c, does not exceed 10% of the 
total level of background for masses up to 2.5 GeV/c 2 . 
Since the accuracy of the JETSET prediction for the 
e + e~ — > 7r + 7r~7T 7r° process is at the 15% level, we con- 
servatively assume that the predictions for non-47r back- 
grounds are accurate to better than 50%. 

The production of a 7r° or a photon with cm. energy 
more than 3 GeV in BB events is kinematically forbid- 
den. We therefore do not expect any background from B 
decays. 

The calculated level of background is about 0.6% in 
the uj mass region and 1.4% in the <f> mass region. The 
3-7T mass spectrum above the <f> for data and background 
events is shown in Fig. 5. The level of background is 
15% at 1.5 GeV/c 2 and 50% at 3 GeV/c 2 . The system- 
atic uncertainty on the background level is about 25% 



below 2 GeV/c 2 . For higher masses the fractional un- 
certainty grows due to an uncertainty in the qq back- 
ground calculation, a model dependence in the simula- 
tion of e + e~ — ► 2-7T7 and e + e~ — > 4-7T7 (the processes 
e + e~ — ► 2tt and e + e~ — > 4-7T have not been measured 
for e + e~ cm. energies above 2 GeV), and the possi- 
ble contribution of unaccounted ISR processes (e + e~ — > 

37T°7, 7r + 7T _ 47T°7, etc.). 

We therefore use two different methods for background 
subtraction. For masses below 1.05 GeV/c 2 , where the 
level of background is low, we subtract the calculated 
background. For higher masses we use the procedure 
of statistical subtraction based on the difference in x 2 
distributions of signal and background events. 

The statistical subtraction procedure is as follows. For 
each mass bin we find the numbers of events with x 2 < 20 
(Ni) and 20 < \ 2 < 40 (N 2 ) and calculate the number 
of signal and background events with \ 2 < 40 as 



N, 



bkg 



(1 _ fi Nl - p Ni 

a- 13 



aN 2 - (1 - a)Ni 
a-P 



(4) 



(5) 



where a,P = Ni/(Ni + N 2 ) for pure signal and back- 
ground events, respectively. Ni and N 2 in Eq. (4-5) do 
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FIG. 8: The 3n mass spectrum for data obtained after the 
statistical background-subtraction procedure. 



FIG. 9: The 3tv mass spectrum for background events ob- 
tained from the statistical background-subtraction procedure 
(points with error bars). The histogram shows the back- 
ground mass distribution expected from MC simulation. 



not contain events e + e 



K + K 7r°7 and e + e 



7T + 7r 7r°7r° processes, which are subtracted from both 
mass distributions bin-by-bin. 

The coefficient (3 is determined from the simulation. 
Its values for the main background processes in four mass 
regions are listed in Table I. It is seen that there is no 
significant dependence of (3 on the mass and that the 
three main background process have consistent values of 
(3. We therefore use the average value (3 = 0.33 ± 0.02 ± 
0.05. The variation of the (3 values for different processes 
was used as an estimate of the systematic error. 

The values of a at the <j) and J/ip masses are extracted 
from data. In the <fi mass region we determine the ra- 
tio Ni/(N\ + N2) for pure signal events from the ex- 
perimental x 2 distribution with subtracted background. 
The value of a at the J/ip mass is measured by another 
method. The numbers of J/ip events for different \ 2 
cuts are determined using a mass-sideband subtraction 
method (see Sec. VIII). The resulting values of a are 
0.879 ± 0.006 ± 0.005 at the <j> mass and 0.882 ± 0.014 at 
the J/ip mass. The systematic error on a at the <f> mass 
is estimated by varying the calculated background level 
by ±25%. For the mass range 1.05-3.00 GeV/c 2 we use a 
linear interpolation between the (p and J/ip values of a. 

The 37r mass spectrum obtained after this statistical 
background subtraction is shown in Fig. 8. The 3-7T mass 
distribution for the background events obtained using 
Eq. (5) is shown in Fig. 9 and compared with the MC sim- 
ulation. The simulation describes well both the shape of 
the background spectrum and the total number of events 
up to at least 2.5 GeV/c 2 . The shapes of the signal and 
background distributions are quite different. 
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FIG. 10: The 37r-mass dependence of the detection efficiency 
obtained from MC simulation. The line is the fit to a second- 
order polynomial. 



V. DETECTION EFFICIENCY 

To first approximation the detection efficiency is de- 
termined from MC simulation as the ratio of the true 37T 
mass distribution computed after and before applying the 
selection criteria. The detection efficiency calculated in 
this way, shown in Fig. 10, is fit to a second-order polyno- 
mial. This efficiency (emc) must be corrected to account 
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TABLE I: The values of p = N{% 2 < 20)/N(x 2 < 40) calculated for different background processes in four mass regions. 
w is the relative contribution of the process to the total background. On the average, each value of j3 is weighted by the 
corresponding value of w. 





1.05 < M 37r < 1.40 


1.4 < M 37T < 2.0 


2.0 < M 3w < 2.5 


2.5 < M 37I < 3.0 


1.05 < M 3n < 3.00 


process 


w(%) 


P 


w(%) 


P 


w(%) 


P 


w(%) 


P 


w(%) 


P 


7r + 7T~7r U 7T U 7 


68 


0.30 ±0.02 


70 


0.33 ±0.02 


67 


0.34 ±0.03 


35 


0.36 ±0.05 


66 


0.32 ±0.01 


7T + 7T~7, /X + ^l~7 


28 


0.33 ±0.10 


20 


0.35 ±0.14 


27 


0.54 ±0.20 


42 


0.29 ±0.19 


26 


0.36 ±0.07 


qq, t + t~ 


4 


0.39 ±0.12 


10 


0.40 ±0.08 


6 


0.50 ±0.20 


23 


0.30 ±0.11 


8 


0.38 ± 0.06 


Average 




0.32 ±0.03 




0.34 ±0.03 




0.40 ±0.07 




0.32 ± 0.09 




0.33 ± 0.02 



for data-MC simulation differences in detector response: 

£ = EMC /n(i + <y, (6) 

where Si are efficiency corrections, for each of several ef- 
fects. These corrections are discussed below and sum- 
marized in Table II. They are determined at the w, <j), 
and J/tjj masses, where the relative level of background 
is small, and a linear interpolation between their values 
is used for the mass ranges between the lo, <p, and J ftp. 

Our preliminary selection contains a cut on the en- 
ergy deposited by charged pions in the calorimeter 
(Eemc/p < 0.9), which is not simulated accurately. The 
momentum dependence of the probability for a pion to 
have Eemc/p > 0.9 is found using events for the process 
e + e~ — > 7T + 7r~7r + 7r~7 selected without cuts on energy 
deposition in the calorimeter. The value of the efficiency 
correction is about 3%. 

The efficiency correction for the background suppres- 
sion cuts is determined from ratios of the number of 
events selected with and without the background sup- 
pression cuts, in data and MC simulation. Wc use events 
with x 2 < 20 from mass ranges near the lo, <f>, and J/4>- 
In data the fraction of signal events rejected by these cuts 
varies from 13% in the u> and </> mass region to 20% at 
J This dependence is reproduced by the simulation. 
The efficiency correction is about 3% for all masses. 

To determine the efficiency correction due to the x 2 < 
40 cut, we fit the 37r mass spectrum in the <fi — oj mass 
range for events with x 2 < 40, 40 < x 2 < 500, and 500 < 
X 2 < 1000 to a sum of simulated signal and background 
distributions with free scale factors. The fit for events 
with 40 < x 2 < 500 is shown in Fig. 11. The data-MC 
simulation discrepancy is estimated by the double ratio 

( , = N slg (A0 < x 2 < 500)/iVMc(40 < X 2 < 500) 
9[ ' A^: 5 (x 2 <40)/iVMc(x 2 <40) 

where Nmc is the number of simulated events and N S i g 
is the number of signal events in data obtained from the 
fit to in the corresponding x 2 interval. We obtain 
#(500) = 1.30 ± 0.04 and, for the 500 < x 2 < 1000 in- 
terval, 5(1000) = 1.29 ± 0.17. For higher x 2 the relative 
signal level is too small to determine its value. From the 
measured values of g(500) and <?(1000) and simulated x 2 
distributions, we calculate the efficiency correction for 
the mass region to be (9 ± 3)%. The quoted error 
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FIG. 11: The 37r mass spectrum for events with 40 < x 2 < 
500. Points with error bars show the data distribution. The 
solid line shows the sum of the distributions for simulated 
e + e~ — > 7r + 7r _ 7r°7 and background processes. The shaded 
histogram shows the distribution for background processes. 

includes uncertainties due to errors on the g values and 
the uncertainty due to events with x 2 > 1000. For these 
events the g value found for the 500 < x 2 < 1000 interval 
is used. In the J/ ijj mass region we use a mass-sideband 
subtraction method to determine the numbers of signal 
events with x 2 < 40 and 40 < x 2 < 500 in data and MC 
simulation. For higher x 2 we do not see a J/i(j signal 
due to large background levels. The obtained value of 
5(500) = 1.12 ± 0.20 agrees with the result for the <fi-Lu 
mass region. Using this number for all x 2 we nn d the 
efficiency correction at the J/ip mass to be (4 ± 6)%. 

The other possible source of data-MC simulation dif- 
ference is 7T° losses due to the merging of electromagnetic 
showers of the two photons from the ir° decay or the loss 
of one of the decay photons. To study the tt° losses we 
perform a kinematic fit for data and simulated events 
with the e + e~ — > ir + ir~ir°~/ hypothesis using the mea- 
sured parameters for only the two charged tracks and 
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FIG. 12: The 37r mass distribution for events selected without 
using the photons from the n° decay, that are reconstructed 
(left) and not reconstructed (right) with our standard kine- 
matic fit procedure. Points with error bars show the data 
distribution. The solid line is a fit result. The shaded his- 
togram is the fitted background contribution. 



the ISR photon. The ir° energy and angles are deter- 
mined as a result of the fit. To suppress background we 
use a very tight cut on the x 2 of the fit and require that 
the total energy of all photons in the event, excluding 
the ISR photon candidate, is greater than 80% of the 
7T° energy found in the fit. The high level of remaining 
background does not allow determination of the efficiency 
correction for masses above the to. Therefore, we restrict 
our study to the w-mass region. The 37r mass spectra 
for selected events reconstructed (x 2 < 10000) and not 
reconstructed with our standard kinematic fit procedure 
are shown in Fig. 12. The fraction of events that are not 
reconstructed is about 30%. In most of these events one 
of the photons from the 7r° decay is outside the polar 
angle range used in our standard selection. The mass 
spectra are fit to a sum of distributions for signal and 
background events. The signal distribution is extracted 
from the simulation. The background spectrum is a sum 
of the simulated distribution for e + e~ — * 7r + 7r _ 7r 7r°7 
events and a second order polynomial with free coeffi- 
cients. The efficiency correction due to ir° losses is de- 
termined to be 5 = eMc/ e ex P — 1 = — (1.9 ± 0.9)%. Here 
e is the fraction of events with a 7r° reconstructed with 
the constrained fit discussed above that, after the stan- 
dard kinematic fit procedure, pass the x 2 < 10000 se- 
lection. To calculate the correction for higher masses we 
must take into account the dependence of the shape of 
the 7T° energy spectrum on the 37r mass. At the u> mass 
we determine the correction as a function of ir° energy 
and convolve this function with the 7r° energy spectra at 
the </> and J/ip masses. The calculated corrections are 
-(1.7 ±0.9)% and -(1.5 ±0.8)% for <j> and J/ip masses, 
respectively. 

We also studied the quality of the simulation of trig- 
ger and background filters used in event reconstruction. 
The corresponding efficiency corrections are listed in Ta- 
ble II. We use the overlap of the samples of events passing 
cither different filters or trigger criteria and the partial 
independence of these filters or triggers to measure their 
efficiency. 
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FIG. 13: Distributions of the ISR photon energy (1st row) and 
polar angle (2nd row) for data (points with error bars) and 
simulation (solid line) events with \ 2 < 20. The left and right 
columns correspond to the mass regions 0.75 < Mz-n < 0.82 
GeV/c 2 and 1.40 < M 37r < 1.80 GeV/c 2 , respectively. Shaded 
histograms show the calculated background contributions. 



The data-MC simulation difference in track losses is 
studied by comparing the ratios of e + e~ — * 7r + 7r~7r + 7r~7 
events with three and four tracks in data and MC simula- 
tion. The difference in data and simulated probabilities 
to lose one of four tracks is found to be (3.6 ± 2.0)%. 
For the case of two tracks we estimate the correspond- 
ing efficiency correction to be (1.8 ± 1.8)%. We increase 
the systematic error to account for the possible depen- 
dence of the correction on track multiplicity and track 
momenta. 

The data-MC simulation difference in the probability 
of photon conversion in the detector material before the 
DCH is studied using e + e~ — ► 77 events and is found to 
be — (0.4 ± 0.2)%. We estimate that the total correction 
for conversion of one of the three photons in an event is 
— (1.0 ± 0.6)%. The fact that part of this correction is 
already included in the correction for 7r° loss is accounted 
for in the determination of this value. 

The event generator for the e + e~ — > 3irj reaction uses 
a model with an intermediate pir state. This model has 
been checked for the lo — <f> mass region in several high 
statistics e + e~ experiments [25-27]. Small deviations 
from the pir model found in the <j> meson decays [27] are 
negligible at our statistical level. Figures 13, 14, and 
15 demonstrate a good agreement between the data and 
simulated distributions over both the kincmatical and dy- 
namical (di-pion invariant masses) parameters for the u> 
mass region. 

In the mass region above the <£> the intermediate state 
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FIG. 14: Distributions of 7r° energy (1st row), charged pion 
momentum (2nd row), and angle between the normal to the 
3ir plane in the three-pion rest frame and the direction of 
the 37T system in the laboratory frame (3rd row) for data 
(points with error bars) and simulation (solid line) events with 
\ 2 < 20. The left and right columns correspond to the mass 
regions 0.75 < M 3n < 0.82 GeV/c 2 and 1.40 < Af 37r < 1.80 
GeV/c 2 , respectively. Shaded histograms show the calculated 
background contributions. 
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lutt becomes noticeable. This additional mechanism was 
studied at the SND experiment [6] for the energy region 
up to 1.4 GeV. It was established that the contribution 
of the lott intermediate state to the total cross section of 
e+e~ — ► 3ir process does not exceed 10%. 

The distributions of different kinematic parameters for 
the mass region from 1.4 to 1.8 GeV/c 2 are shown in 
Figs. 13 and 14. For these parameters the data distri- 
butions agree with the simulated ones. In the tt + tt~ 
mass spectra (Fig. 15) for data, a narrow peak near the 
uj mass is seen. The fraction of events in this peak is 
about 10% for 1.1 < M 3w < 1.4 GeV/c 2 and 6% for 
1.4 < < 1.8 GeV/c 2 in agreement with measure- 

ments of SND [6]. We calculate the detection efficiency 
for simulated events with 0.77 < M n + n - < 0.80 GeV/c 2 



FIG. 15: Distributions of ty + tt~ (left column) and^ir (right 
column) invariant masses for data (points with error bars) 
and simulation (solid line) for events with \ 2 < 20. The 
rows correspond to the mass regions: 0.75 < Mz-k < 0.82 
GeV/c 2 (1st row), 1.00 < M [i7V < 1.04 GeV/c 2 (2nd row), 
1.10 < M 3n < 1.40 GeV/c 2 (3rd row), and 1.40 < M Sn < 1.80 
GeV/c 2 (4th row). Shaded histograms show the calculated 
background contributions. 



and find its difference with detection efficiency for the 
full M w +„- mass range to be (1.6 ± 3.9)%. 

There are also noticeable deviations from the MC sim- 
ulation in the distribution of the tt^tt invariant mass for 
the mass region 1.4 < il/ 37r < 1.8 GeV/c 2 (Fig. 15), which 
are manifest in the shift of the visible p-mcson mass po- 
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sition and a bump at a mass above 1 GcV/c 2 . A possible 
explanation of this difference is the appearance at higher 
masses of the uj" — > p'n transition, which interferes with 
the pir amplitude. We study the dependence of the detec- 
tion efficiency on tt + tt° mass and find it to be constant 
in the mass range 0.65 < M„+„o < 1.25 GeV/c 2 . Wc 
conclude that the use of the pit model for the simulation 
of e + e~ — > 3-7T7 does not lead to any significant errors in 
the determination of the detection efficiency. 

Efficiency corrections Si are summarized in Table II. 
The total efficiency corrections S calculated from 1 + 5 = 
11(1 + 5i) are listed in the last row of the table. The 
correction is 14% near the ui and <f>, and 11% at the J/ip 
mass. 



where 03^(771) is the Born cross section for e + e~ — > 3ir, 
dL/dm is the so-called ISR differential luminosity, e is 
the detection efficiency as a function of mass, and R is a 
radiative correction factor accounting for the Born mass 
spectrum distortion due to emission of several photons 
by the initial electron and positron. The ISR luminos- 
ity is calculated using the total integrated luminosity L 
and the probability density function for ISR photon emis- 
sion (Eq. (2)): 

di a ( . „ , 1 + C J\ 2m . . 

— = — 2 - 2x + x 2 ) log — — - x 2 C — L. (8) 
dm ttx \ 1 — G / s 



VI. 



DISTRIBUTION 



In order to determine the peak cross sections for e + e~ 
annihilation into u> and <f> mesons and the resonance 
parameters of excited u> states, we fit the background- 
subtracted 37r invariant-mass spectrum. The mass spec- 
trum is described by the following function: 



dN 
dm 



^{m)— Re, 
dm 



(7) 



Here x = 1 — m 2 /s, y/s is the e + e~ cm. energy, C = 
cos9 , and 9 determines the range of polar angles in the 
cm. frame: O < 6 1 < 180° - 6 for the ISR photon. 
In our case 9q is equal to 20°, since we determine the 
detector efficiency using the simulation with 20° < 6* 7 < 
160°. 



37r can be written 
as the sum of the contributions of four resonances: 



12-7T . . ^— v 

cr37r(m) = -^F p7T {m) ) j 



T v m^ 2 ^B{y 



e+e~)B{V 3tt) e l<t > v 



V=u 



D v {m) 



\fFp~Jmv) 



(9) 



where my and IV are the mass and width of the reso- 
nance V, 4>v is its phase, B(V — > e + e _ ) and B(V — > 37r) 
are the branching fractions of V into e + e~ and 37r, 

Dy(m) = my — m 2 — imTy(rn), Yy(rri) = Tf(m). 

f 

Here Tf(m) is the mass-dependent partial width of the 
resonance decay into the final state /, and IV (my) = 
TvB(V — > /). The mass-dependent width for the u> and 
cj) mesons has been calculated taking into account all sig- 
nificant decay modes. The corresponding formulae can 
be found, for example, in Ref. [25]. We assume that 
V — > 37r decay proceeds via the pn intermediate state, 
and Fpn (m) is the 3-7T phase space volume calculated un- 
der this hypothesis. The formula for F, m calculation can 
be found in Ref. [25]. 

The radiative correction factor was determined using 
Monte Carlo simulation (at the generator level, with 
no detector simulation). The 37T mass spectrum was 
generated both using only the pure Born amplitude 
of the e + e~ — > ir+ir^ir ^ process and using a model 



with higher-order radiative corrections included with the 
structure function method. With the cut on the invari- 
ant mass of the 3tt^/ system, M 37r7 > 8 GeV/c 2 , used in 
our simulation, no significant difference is found between 
these two spectra. Therefore the radiative correction fac- 
tor is evaluated as the ratio of the total cross section with 
M 37r7 > 8 GeV/c 2 to the Born cross section and is found 
to be close to unity, R = 0.9994. The theoretical un- 
certainty in the radiative correction calculation with the 
structure function method does not exceed 1% [19]. The 
radiative correction factor does not include the correc- 
tions due to leptonic and hadronic vacuum polarization. 
Here we follow the generally accepted practice [28] of 
including the vacuum polarization correction in the reso- 
nance electronic width. The probability density function 
(PDF) for the 37r mass spectrum as expressed in Eq. (7) 
needs to be convolved with the detector resolution func- 
tion in order to fully characterize the experimental mass 
distribution found in the data. The detector resolution 
function is obtained using MC simulation of the detec- 
tor response. Figure 16 shows the distribution of the 
difference between measured and true 3-7T mass for sim- 
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TABLE II: The values of different efficiency corrections 5, for uj, (f>, and J/ijj mass regions. 



effect 


5i(m u ),% 


Si(m#),% 




Eemc/p < 0.9 cuts 


+2.9 ±0.3 


+2.9 + 0.3 


+2.5 + 0.3 


nfi^KormiTin rpippfinn r"iii"Q 


+3.0 ± 0.4 


+2.4 ± 0.7 


+2.9 ±1.9 


X 2 < 40 cut 


+9 + 3 


+9 + 3 


+4 + 6 


7r° loss 


-1.9 + 0.9 


-1.7 + 0.9 


-1.5 ±0.8 


trigger and filters 


+0.0 + 0.4 


+0.1 + 1.0 


±2.2 ±2.0 


track loss 


+1.8+1.8 


+1.8 + 1.8 


+1.8 ± 1.8 


photon conversion 


-1.0 + 0.6 


-1.0 + 0.6 


-1.0 + 0.6 


total 


+14 + 4 


+14 + 4 


±11 ±7 




FIG. 16: The distribution of the difference between measured 
and true 37r mass for simulated events with \ 2 < 40 from the 
mass region near the u. The curve is a fit to a triple-Gaussian 
function. 



ulated events with x 2 < 40 in the mass region near the 
uj. For each Sir mass region, the distribution is fit with 
a triple-Gaussian function. The resolution depends on 
the 3-7T mass and is about 6, 7, and 9 MeV/c 2 at the uj, 
<f> and J/ip masses, respectively. The determination of 
the resolution function is performed in the u>, <f> and J /ip 
mass regions, where the available Monte Carlo simulation 
statistics are high enough. For the mass region between 
the uj and <f> we use a linear interpolation. Since no nar- 
row peaks are present in the mass region between the <f> 
and the J/ip, it is not critical to know a very detailed 
resolution function, and therefore an average resolution 
function is used over the mass range between 1.05 and 
3.00 GcV. 

A binned maximum likelihood fit is used to fit the 
3tt mass spectrum in data. The bin width is chosen 
to be 2.5 MeV/c 2 for the uj — cf> mass range and 25 
MeV/c 2 for masses above the c\>. The free parameters 
in the fit are the products of the branching fractions 



^ 300 

o 

> 
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FIG. 17: The background-subtracted 37r mass spectrum for 
masses between 0.70 and 1.05 MeV/c 2 (upper plot) and for 
masses from 1.05 to 1.80 MeV/c 2 (lower plot). The curves 
are the result of the fit described in the text. 
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B(V — > e + e~)B(V — > 371"), the masses my for all four 
resonances (V = u,cf),Lv',uj"), and the widths for lo' and 
lo" . The width of the u> meson is fixed at the value of 
(8.68±0.13) MeV/c 2 obtained recently in the CMD2 [29] 
and SND [25] experiments. The width of the cf> meson 
is fixed at the PDG value. The relative phase between 
the lo and <fi amplitudes, (f)^ = (163 ± 7)°, is taken from 
Ref. [25]. The phases of w, ui', and lo" are fixed at values 
of 0°, 180°, and 0° [30]. Our fitting function docs not take 
into account the contribution of the e + e~ — > lot:^ — > 37r7 
process which proceeds via excited p states. In Ref. [6] 
it is shown that this mechanism does not change signifi- 
cantly the parameters for the to' and lo" resonances. The 
fitted mass region is restricted to masses below 1.8 GeV. 
To be described properly, data in the higher mass range 
would require a more complicated function, which would 
take into account both the resonant and the non-resonant 
37T production. There are no reliable models available in 
the literature, and therefore our results on the parame- 
ters of the lo 1 and lo" states, obtained in this fit, should 
only be considered a first approximation. 

In order to account for a possible resolution differ- 
ence between data and simulation, the resolution func- 
tion determined from simulation is modified by adding or 
subtracting quadratically an additional og to all sigmas 
of the triple-Gaussian function. Technically, a squared 
sigma <Jq is used as a free parameter (with negative val- 
ues allowed). Two a G parameters are used, one for lo and 
another for <p and higher masses. 

The fit result is shown along with the data in Fig. 17. 
The resulting parameters obtained from the fit (x 2 /dof = 
146/148) are the following: 

B{lo -> e+e-)B{Lo -> 3tt) = (6.70 ± 0.06 ± 0.27) x 10~ 5 , 
m u - toP° g = -(0.2 ± 0.1) MeV/c 2 , 



a Guj = (0.9 ± 1.6) MeV 2 /c 4 



B( 



e+e-)B(<j) -► 3tt) = (4.30 ± 0.08 ± 0.21) x 10 
m 4> - m l BG = -(°- 6 ± °-2) MeV/c 2 , 



-5 



Gtj> 



-(3.2 ±2.6) MeV 2 /c 4 , 
B(lo' -> e+e-)B(u' -> 3tt) = (0.82 ± 0.05 ± 0.06) x 10~ 6 , 
M u > = (1350 ± 20 ± 20) MeV/c 2 , 
T u > = (450 ± 70 ± 70) MeV/c 2 , 
B{lo" -> e+e-)B(uj" -> 3tt) = (1.3 ± 0.1 ± 0.1) x 10~ 6 , 
M u „ = (1660 ± 10 ± 2) MeV/c 2 , 
Y u „ = (230 ± 30 ± 20) MeV/c 2 . 

The quoted errors correspond to the statistical and sys- 
tematic uncertainties, respectively. The systematic error 
for B(V — > e + e~)B(V — > Sir) includes a statistical error 
from simulation, the error on the efficiency correction 
(Table II), 1.2% uncertainty in the luminosity, 1% the- 
oretical uncertainty on the radiative correction, a back- 
ground subtraction uncertainty (0.4% at lo and 0.6% at 
(f>), and an uncertainty arising from errors on r w , F^, and 
(1% at lo and 2.8% at (f). The systematic errors on 
the masses and widths of the lo' and lo" mesons are due 



to the background-subtraction uncertainty and the errors 
on T u , r<£, and <j>^. 

The fitted values B(V -► e+ e -)B(V -> 3tt) for the w 
and <fi mesons are in reasonable agreement with the corre- 
sponding world average values [5], B(lo — > e + e~)B(to — > 
3tt) = (6.35 ± 0.10) x 10~ 5 and B(<j) -> e+e~)B{(j) -> 
3tt) = (4.52 ±0.19) x 10~ 5 . The observed peak positions 
of both the lo and 4> are shifted to lower masses relative to 
their PDG values. The shifts are about (0.3-0.6) xl0~ 3 
of the mass values. The fitted values of a G parameters 
have large statistical uncertainties, and lead to a change 
in the simulated resolution of (2-3)%. 

The fitted masses and widths of the lo' and lo" mesons 
can be compared with the estimates of these parameters 
by the PDG [5]: M u > = 1400 - 1450 MeV/c 2 , IV = 
180 - 250 MeV/c 2 , M u » = 1670 ± 30 MeV/c 2 , = 
315 ± 35 MeV/c 2 . The PDG data are based on small 



data samples for e H 



— > lo' ,lo" — > 37T, lotttt [6, 7, 31], 
pp — > lo'it — > o;7r 7r 7r [32], and ir~p — > Lo"n — > Lorjn [33] 
reactions. We present a new measurement of the lo' and 
lo" parameters based on a significantly larger data sam- 
ple for the e + e~ — > lo',lo" — > 3ir reaction. From the 
measured values of B(V — * e + e~)B(V — ► 3ir), the elec- 
tronic widths of lo' and lo" can be estimated. Assuming 
that B{lo' — > 37r) w 1 and B{lo" — > 37r) w 0.5 we derive 
that T(uj' e+e~) w 370 eV and r(w" -> e+e^) w 570 
cV. The large values of these widths, comparable with 
T(lo — * e + e _ ) w 600 eV, are in disagreement with expec- 
tations of the quark model, which predicts at least one 
order of magnitude lower values for the electronic widths 
for the excited meson states (see, for example, Ref. [34]). 



VII. MEASUREMENT OF THE e+e" 
CROSS SECTION 

The cross section for e + e~ — > 7r + 7r~7r°, in the energy 
(Vs') range between 1.05 and 3.00 GeV, is calculated 
from the 3tt mass spectrum using 



<T37r(m) 



(dJV/dm) et 
sRdL/dr 



(10) 



where m = ^/s* is the 3n invariant mass, and 
(dN I dm) corr is the mass spectrum corrected for reso- 
lution effects. 

The resolution-corrected mass spectrum is obtained by 
first subtracting the events with an actual 37r invariant 
mass outside the 1.05-3.00 GcV/c 2 region (tails of the 
cf> and J/tp mass distribution). The number of 0- meson 
events with measured mass above 1.05 GeV/c 2 is esti- 
mated from simulation. We subtract 10 ± 5 events from 
the first mass bin (1.05-1.075 GeV/c 2 ). The number of 
J/ip events contributing to the mass region under study 
is found to be 1 ± 1. Second, the detector resolution is 
deconvolved by using a migration matrix A that gives 
the probability that an event with true mass in bin j is 
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FIG. 18: 
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7r + 7r 7T° cross section measured in this work (filled circles), by SND (open circles), and DM2 (open 



actually reconstructed in bin i: 



(11) 



The inverse of this migration matrix {A^ 1 ) is then ap- 
plied to the measured spectrum. In our case where the 
spectrum has no narrow structures and the smearing is 
small, we can neglect the mass dependence of the reso- 
lution and use the average resolution function to build 
the migration matrix. This allows a determination of the 
inverse matrix that is robust against statistical fluctu- 
ations. The resolution function is obtained from simu- 
lation and takes into account our background subtrac- 
tion procedure. In practice, the only elements of the mi- 
gration matrix that are significantly different from zero 
are the diagonal elements, which are around 0.84, and 
the elements next to the diagonal, which are around 
0.08. The resolution correction procedure changes in- 
significantly the shape of the mass distribution but leads 
to an increase in the errors (by «20%) and their corre- 
lation. The number of events in each bin changes by no 
more than half its statistical uncertainty. In particular, 
near 37r mass of 1.6 GeV this change is less than 3%. The 
number of events for each mass bin is listed in Table III. 
The quoted errors correspond to the statistical (includ- 
ing the contribution of the background subtraction for 
the e + e~ — > K + K~tt°j and e + e~ — ► 7r + 7r~7r°7r° pro- 
cesses) and systematic (due to the a and (3 parameters) 
uncertainties. 

Table III also contains the values of the detection effi- 



ciency and the ISR differential luminosity calculated ac- 
cording to Eq. (8). The errors on the detection efficiency 
include the uncertainties on the efficiency correction and 
the statistical errors from the simulation. 

The calculated cross section is shown in Fig. 18 and 
listed in Table III. The quoted errors are statistical and 
systematic. The latter includes the systematic error con- 
tributions from the number of events, and uncertainties 
in the detection efficiencies and in the calculation of the 
radiative correction. Note that the systematic errors for 
different mass bins are fully correlated. 

Most of experiments at low energy measure so-called 
"dressed" cross sections (see, for example, Ref. [6]), which 
include the vacuum polarization corrections. Since our 
radiative correction factor R does not take into account 
vacuum polarization we also measure the "dressed" cross 
section. A comparison of our measurements with other 
e + e~ data is shown in Fig. 18. Our values are in good 
agreement with SND measurements, but significantly ex- 
ceed DM2 results. 



VIII. MEASUREMENT OF THE J/ip 
BRANCHING FRACTION 



3tt 



The 3-7T mass spectrum for selected events in the J/ip 
mass region is shown in Fig. 19. The small width of the 
J/tp resonance leads to negligible peaking background. 
For example, e + e~ — > J/tjj^f — > K + K ~ir a j events recon- 
structed under the 37T7 hypothesis have a Sir invariant 
mass in the range 2.8 to 3.0 GcV/c 2 . Since the nonreso- 
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TABLE III: e + e~ cm. energy (VF), number of selected events after 37r-mass-resolution correction (N corr ), detection efficiency 
(e), differential ISR luminosity (L), and measured cross section (a) for e + e~ — > tv + tv~tv {) . All values are calculated for 25 MeV 
bin size. The quoted errors are statistical and systematic. 



V s' 
(GeV) 


N corr 


£ 
(%) 


L 

(nb ) 








(nb) 








Vs 7 
(GeV) 


N corr 


e 
(%) 


L 

(nb ) 


a 
(nb) 


1.0625 


40 ± 13 ± 3 


9.8 ± 0.4 


254 


1 


61 


± 





.53 


± 





,13 


2 


,0375 


22 ± 9 ± 2 


10.2 ± 0.5 


511 


0.41 ± 0.17 ± 0.04 


1.0875 


94 ± 16 ± 4 


9.8 ± 0.4 


260 


3. 


68 


± 





.61 


± 





.21 


2 


,0625 


23 ± 8 ± 1 


10.2 ± 0.5 


518 


0.43 ± 0.15 ± 0.02 


1.1125 


114 ± 16 ± 3 


9.8 ± 0.4 


266 


4 


35 


± 





,63 


± 





.21 


2 


,0875 


33 ± 9 ± 1 


10.2 ± 0.5 


525 


0.61 ±0.17 ±0.04 


1.1375 


133 ± 18 ± 4 


9.9 ± 0.4 


272 


4 


96 


± 





,66 


± 





,25 


2 


,1125 


28 ± 9 ± 1 


10.2 ± 0.5 


532 


0.52 ± 0.16 ± 0.03 


1.1625 


115 ± 17 ± 4 


9.9 ± 0.4 


278 


4 


17 


± 





,62 


± 





22 


2 


1375 


24 ± 8 ± 1 


10.2 ± 0.5 


539 


0.44 ± 0.15 ± 0.03 


1.1875 


150 ± 18 ± 4 


9.9 ± 0.4 


285 


5. 


31 


± 





,65 


± 





.25 


2 


,1625 


30 ± 9 ± 1 


10.2 ± 0.5 


547 


0.54 ± 0.16 ± 0.03 


1.2125 


150 ± 19 ± 4 


9.9 ± 0.4 


291 


5. 


19 


± 





,65 


± 





.26 


2 


,1875 


24 ± 8 ± 1 


10.2 ± 0.5 


554 


0.43 ± 0.15 ± 0.03 


1.2375 


198 ± 21 ± 5 


9.9 ± 0.4 


297 


6. 


69 


± 





,71 


± 





.31 


2 


,2125 


20 ± 8 ± 1 


10.2 ± 0.5 


561 


0.36 ± 0.14 ± 0.02 


1.2625 


185 ± 20 ± 4 


10.0 ± 0.4 


304 


6. 


11 


± 





,65 


± 





.27 


2 


,2375 


15 ± 8 ± 1 


10.2 ± 0.5 


569 


0.26 ± 0.13 ± 0.02 


1.2875 


151 ± 19 ± 5 


10.0 ± 0.4 


310 


4 


88 


± 





,62 


± 





.25 


2 


,2625 


23 ± 8 ± 1 


10.2 ± 0.5 


576 


0.39 ± 0.14 ± 0.03 


1.3125 


214 ± 21 ± 4 


10.0 ± 0.4 


316 


6. 


75 


± 





,65 


± 





.31 


2 


,2875 


21 ± 8 ± 1 


10.2 ± 0.5 


584 


0.36 ± 0.14 ± 0.02 


1.3375 


115 ± 18 ± 6 


10.0 ± 0.4 


323 


3. 


54 


± 





,55 


± 





.23 


2 


,3125 


23 ± 8 ± 1 


10.2 ± 0.5 


591 


0.39 ± 0.14 ± 0.03 


1.3625 


160 ± 19 ± 4 


10.0 ± 0.4 


329 


4 


85 


± 





,56 


± 





.24 


2 


,3375 


5 ± 7± 1 


10.1 ± 0.5 


599 


0.09 ± 0.11 ± 0.01 


1.3875 


136 ± 18 ± 4 


10.1 ± 0.4 


335 


4 


01 


± 





,53 


± 





21 


2 


,3625 


28 ± 9 ± 1 


10.1 ± 0.6 


606 


0.46 ± 0.14 ± 0.03 


1.4125 


150 ± 19 ± 5 


10.1 ± 0.4 


342 


4 


36 


± 





,54 


± 





.24 


2 


,3875 


9 ± 7± 1 


10.1 ± 0.6 


614 


0.15 ± 0.11 ± 0.01 


1.4375 


116 ± 17 ± 5 


10.1 ± 0.4 


348 


3. 


31 


± 





,49 


± 





,20 


2 


,4125 


23 ± 8 ± 1 


10.1 ± 0.6 


621 


0.37 ± 0.12 ± 0.02 


1.4625 


156 ± 19 ± 4 


10.1 ± 0.4 


355 


4 


36 


± 





,54 


± 





21 


2 


,4375 


20 ± 8 ± 1 


10.1 ± 0.6 


629 


0.32 ± 0.12 ± 0.02 


1.4875 


144 ± 19 ± 3 


10.1 ± 0.4 


361 


3. 


93 


± 





,52 


± 





.20 


2 


,4625 


9 ± 6 ± 1 


10.1 ± 0.6 


637 


0.14 ± 0.10 ± 0.01 


1.5125 


172 ± 20 ± 4 


10.1 ± 0.4 


368 


4 


60 


± 





,53 


± 





.22 


2 


,4875 


17 ± 7 ± 1 


10.0 ± 0.6 


645 


0.26 ± 0.11 ± 0.02 


1.5375 


146 ± 19 ± 3 


10.2 ± 0.4 


374 


3. 


83 


± 





,50 


± 





19 


2 


,5125 


16 ± 7 ± 1 


10.0 ± 0.6 


653 


0.24 ± 0.11 ± 0.02 


1.5625 


220 ± 22 ± 5 


10.2 ± 0.4 


381 


5. 


69 


± 





,56 


± 





,28 


2 


,5375 


7 ± 6 ± 1 


10.0 ± 0.6 


660 


0.11 ± 0.09 ± 0.01 


1.5875 


213 ± 22 ± 5 


10.2 ± 0.4 


387 


5. 


39 


± 





,56 


± 





.27 


2 


,5625 


16 ± 7 ± 1 


10.0 ± 0.6 


668 


0.24 ± 0.10 ± 0.02 


1.6125 


250 ± 23 ± 5 


10.2 ± 0.5 


394 


6. 


24 


± 





,58 


± 





.31 


2 


,5875 


3 ± 6 ± 1 


10.0 ± 0.6 


676 


0.05 ± 0.08 ± 0.01 


1.6375 


218 ± 22 ± 5 


10.2 ± 0.5 


401 


5. 


33 


± 





,55 


± 





.27 


2 


,6125 


21 ± 7 ± 1 


9.9 ± 0.6 


684 


0.30 ± 0.10 ± 0.02 


1.6625 


201 ± 21 ± 4 


10.2 ± 0.5 


407 


4 


84 


± 





,50 


± 





.24 


2 


,6375 


6±6± 1 


9.9 ± 0.6 


693 


0.08 ± 0.08 ± 0.01 


1.6875 


107 ± 17 ± 3 


10.2 ± 0.5 


414 


2. 


52 


± 





,40 


± 





.14 


2 


,6625 


13 ± 7 ± 1 


9.9 ± 0.6 


701 


0.18 ± 0.10 ± 0.01 


1.7125 


100 ± 16 ± 2 


10.2 ± 0.5 


421 


2. 


32 


± 





,36 


± 





12 


2 


,6875 


5 ± 7 ± 1 


9.9 ± 0.6 


709 


0.07 ± 0.09 ± 0.01 


1.7375 


62 ± 14 ± 2 


10.2 ± 0.5 


427 
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42 
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,31 


± 





,08 


2 


,7125 


9 ± 6 ± 1 


9.8 ± 0.6 


717 


0.12 ± 0.09 ± 0.01 


1.7625 


72 ± 14 ± 2 


10.2 ± 0.5 


434 


1 


63 


± 





,31 


± 





.09 


2 


,7375 


2 ± 5 ± 1 


9.8 ± 0.6 


725 


0.03 ± 0.08 ± 0.01 


1.7875 


62 ± 13 ± 2 


10.2 ± 0.5 


441 


1 


38 


± 





,29 


± 





.07 


2 


,7625 


-1 ± 5 ± 1 


9.8 ± 0.6 


734 


-0.01 ± 0.07 ± 0.01 


1.8125 


52 ± 12 ± 1 


10.2 ±0.5 


448 


1 


14 


± 





,26 


± 





.06 


2 


,7875 


0±5± 1 


9.8 ±0.6 


742 


0.00 ±0.07 ±0.01 


1.8375 


46 ± 13 ± 3 


10.2 ±0.5 


455 





99 


± 





,27 


± 





.07 


2 


,8125 


13±7±1 


9.7 ±0.6 


751 


0.17 ±0.09 ±0.01 


1.8625 


64 ± 12 ±2 


10.3 ±0.5 


462 


1 


35 


± 





,26 


± 





.07 


2 


,8375 


4± 6± 1 


9.7 ±0.6 


759 


0.05 ±0.08 ±0.01 


1.8875 


50 ± 12 ± 2 


10.3 ±0.5 


468 


1 


04 


± 





,24 


± 





.07 


2 


,8625 


6±5± 1 


9.7 ±0.6 


768 


0.08 ±0.07 ±0.01 


1.9125 


39 ± 11 ±2 


10.3 ±0.5 


475 





80 


± 





,22 


± 





.06 


2 


,8875 


6±5± 1 


9.6 ±0.6 


776 


0.08 ±0.07 ±0.01 


1.9375 


41±10±1 


10.3 ±0.5 


482 
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± 





,21 


± 





.05 


2 


,9125 


6±5± 1 


9.6 ±0.6 


785 


0.08 ±0.07 ±0.01 


1.9625 


32 ± 9 ± 1 


10.3 ±0.5 


489 





64 


± 





19 


± 





.04 


2 


,9375 


20 ± 7 ± 1 


9.6 ±0.6 


794 


0.26 ± 0.09 ± 0.02 


1.9875 


15±8±1 


10.3 ±0.5 


496 





30 


± 





16 


± 





.03 


2 


,9625 


-1±5± 1 


9.5 ±0.6 


802 


-0.01 ±0.06 ±0.01 


2.0125 


31 ±9± 1 


10.2 ±0.5 


503 





60 


± 





18 


± 





.04 


2 


,9875 


16±6±1 


9.5 ±0.6 


811 


0.20 ±0.08 ±0.01 



signal region (3.0 < M 3w < 3.2 GeV/c 2 ) and in the side- 
bands (2.9 < M 37r < 3.0 GeV/c 2 and 3.2 < M 3w < 3.3 
GeV/c 2 ). 

The Monte Carlo simulation of the number of events 
in the signal and sideband regions is used to estimate a 
detection efficiency of emc = 0.101±0.002. The data- 
MC simulation differences discussed earlier are used to 
correct the former efficiency value by (11 ± 7)%. 

The simulation uses the pn model of J/tp — > 37r decay. 
In order to check the model dependence of the detec- 
tion efficiency, the Dalitz plot for events in the J/ ip peak 
(Fig. 20) is analyzed. It is seen that the main mechanism 
for J/iJj — > 37r decay is pir. There is, however, a differ- 
ence between the data and simulated plots (an absence 



TABLE IV: N 3 i gna i and N s ide are the numbers of selected 
events in the signal region (3.0 < < 3.2 GeV/c 2 ) and the 
sidebands (2.9 < M 37r < 3.0 and 3.2 < M 3lr < 3.3 GeV/c 2 ), 
respectively. 





N signal 


Nside 


signal -^V side 


data 


1023 


103 


920±34 


MC 


1825 


13 


1812±43 



nant background is small and well described by a linear 
function, a mass-sideband subtraction method is used to 
determine the number of J/ip events. Table VIII shows 
the numbers of data and simulated 3777 events in the 
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following product can be determined: 
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FIG. 19: 3ir mass spectrum for selected e + e — > J/4>~f 
7r + 7r~7r°7 events. 
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FIG. 20: The Dalitz plot for J/ip — > 37r candidates in data 
(left) and simulation (right). 



T(J/ip -► e + e-)B{J/ip -► 3tt) 
= (0.122 ± 0.005 ± 0.008) keV. 

The systematic error includes the uncertainties on the 
detection efficiency, the integrated luminosity, and the 
radiative correction. 

The most precise measurement of the electronic width 
was made in the analysis of e + e~ — ► J/tp"f — > /i + /i~7 
by BABAR [8]: T(J/ip -> e+e") = (5.61 ± 0.20) keV. 
Using the latter measurement, the J/ip — > 37r branching 
fraction is calculated to be 

B{J/ip -► 3tt) = (2.18 ± 0.19)%, 

which is in substantial disagreement (~ 3er) with the 
world average value (see Sec. I) of (1.47 ± 0.13)%, but 
agrees with the result from the BES collaboration [16]: 
B(J/ip -» 3tt) = (2.10 ± 0.12)%. 



IX. SUMMARY 

The process e + e~ — > 7r + 7r _ 7r°7 was studied for the 3-7T 
invariant masses up to 3 GeV/c 2 and at the J/ip mass. 
From the measured 3tt mass spectrum we obtained the 
e + e~ — > 7T + 7r _ 7r° cross section for the 1.05 < \fs' < 3 
GeV energy range. The results are in agreement with 
the SND measurement [26] for \fs' < 1.4 and signifi- 
cantly exceed the DM2 data [7] in the 1.4 < y/s 7 < 2.2 
range. The e + e~ — ► ir + ir~ir° cross section in the energy 
range up to 1.8 GeV is described well by a sum of the 
contributions of four isoscalar resonances: w, <p, a/, and 
cu" . From the fit of the 37r mass spectrum we obtained 
the following parameters for these resonances: 



of events in the center of the Dalitz plot for data) , which 
can be a manifestation of negative interference with the 
contribution of intermediate states other than the pir [35] . 
The influence of this difference on the detection efficiency 
is studied by excluding events located in the center of the 
Dalitz plot in the simulated sample, and recomputing the 
detection efficiency. The result is a (1.1±0.6)% rise in ef- 
ficiency. This correction is included with a systematic 
error of 1.1% in the final calculation of the detection ef- 
ficiency, which is determined to be 0.092 ± 0.006. 

The cross section for e + e~ — > J/ip~f — > 7r + 7r~7r°7 for 
20° <6~< 160° is calculated as 



ct(20° < 6 1 < 160°) 



signal side 

eRL 



= (112±4±8) fb. 



The radiative-correction factor R = a jcJBorni of 1.005 ± 
0.002±0.010 used here, is obtained from a MC simulation 
at the generator level (no detector simulation). The total 
integrated luminosity for the data sample is (89.3 ± 1.1) 
fb -1 . From the measured cross section and Eq. (3), the 



B(cu 
B{4>- 
B{lo' - 



B(u>" 



e~)B{u -> 3tt) = (6.70 ± 0.06 ± 0.27) x 10~ 5 , 
e-)B(4> -> 3tt) = (4.30 ± 0.08 ± 0.21) x 10~ 5 , 
■e-)B(w' -> Sir) = (0.82 ± 0.05 ± 0.06) x lO" 6 , 
M u > = (1350 ± 20 ± 20) MeV/c 2 , 

Fy = (450 ± 70 ± 70) MeV/c 2 , 
3+e-)B(uj" -> 3tt) = (1.3 ± 0.1 ± 0.1) x 10~ 6 , 

M u » = (1660 ± 10 ± 2) MeV/c 2 , 

T u » = (230 ± 30 ± 20) MeV/c 2 . 



The electronic widths of u>' and lu" corresponding to 
these resonance parameters, T(oj' — > e + e~) w 370 eV 
and T(lli" — > e + e~) r* 570 cV, are comparable with the 
cj(782) electronic width, in disagreement with expecta- 
tions of the quark model (see, for example, Ref. [34]). 

From the measured number of events in the e + e~ — > 
J/ipj — ► 7r + 7r^7r°7 reaction we determine 



T(j/ip 



e e 



~)B(J/ip -> 3tt) 



(0.122 ± 0.005 ±0.008) keV. 
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Dividing this value by T(J/ip -»■ e+e") = (5.61 ± 0.20) 
keV [8] we obtain B(J/tp -> 3tt) = (2.18 ± 0.19)%, which 
is in ~*3cr disagreement with the world average value of 
(1.47 ± 0.13)% (see Sec. I), but agrees with the recent 
result from the BES Collaboration [16]: (2.10 ± 0.12)%. 
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